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Using 226 million T(4S) -» BB events collected with the BABAR detector at the PEP-II e + e" 
storage ring at the Stanford Linear Accelerator Center, we measure the branching fraction for 
B° -> D°K + ir-, excluding B° -> D*~ K+ , to be B{B° -> D°K + ^-) = (88 ± 15 ± 9) x 10~ 6 . 
We observe B° -> D°K* (892)° and 5° -> D|(2460)~ K + contributions. The ratio of branching 
fractions B(B° -> D*-K+)/B(B° D*~-k + ) = (7.76 ± 0.34 ± 0.29)% is measured separately. The 
branching fraction for the suppressed mode B° -> D°K + tt' is B(B° -> D°K + iv~) < 19 x 10~ 6 at 
the 90% confidence level. 



PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 

A theoretically clean method for measuring the an- 
gle 7 = &vg{-V ud V* b /V cd y* b ) in the unitarity triangle of 
the Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing 
matrix [1] in the Standard Model of particle physics uti- 
lizes decay modes of the type B — > DK. Several meth- 
ods have been proposed [2-4] to extract 7 from these 
decays using interference effects between b — > ucs and 
b — > cus processes. However, the b — > ucs amplitude is 
suppressed by a color factor in addition to the CKM fac- 
tor IV^bV^/V^&V^J ~ 0.4, and the extraction of 7 with 
methods in Rcf. [2, 3] is subject to an eight-fold ambigu- 
ity due to unknown strong phases. 

Three-body B — ► DKtt decays have been proposed [5, 
6] as an alternative method for measuring 7. In these 
modes, the CKM-suppressed b — ► ucs processes include 
color-allowed diagrams; thus larger decay rates and more 
significant CP violation effects are possible. In addition, 
a DKtt Dalitz plot analysis can resolve the strong phase 
and reduce the ambiguity to two-fold, similar to Rcf. [4] . 
The sensitivity to 7 in these decays is determined by the 
size of the overlapping b — ► cus and b — > ucs amplitudes 
in the Dalitz plot. 

In this Letter, we report the measurements of 
the branching fraction for the CKM-favored B° — > 
D°K + n~ [7] decay and dominant resonance contribu- 
tions, and the search for the CKM-suppressed B° — > 
D°K + tt~ decays. The flavor of the B meson is tagged 
by the charge of the prompt kaon. The favored mode 
has been previously observed through its dominant res- 
onances D*~K+ [8] and D°K*(892)° [9]. Since D*~ K+ 
occupies only a very small region of the allowed phase 
space, we treat it separately and measure the ratio r = 
B(B° -> D*-K+)/B(B° -> D*-n+), which can be used 
to test factorization and flavor-SU(3) symmetry. 

Signal events are selected from 226 million BB pairs 
collected with the BABAR detector [10] at the PEP-II 
asymmetric-energy storage ring. Charged tracks are de- 
tected by a five-layer silicon vertex tracker and a 40-laycr 
drift chamber. Hadrons are identified based on the ion- 



ization energy loss in the tracking system and the opening 
angle of the Cherenkov radiation in a ring-image detec- 
tor [11]. Photons are measured by an electromagnetic 
calorimeter. These systems are mounted inside a 1.5-T 
solenoidal super-conducting magnet. 

The D° candidate is reconstructed through K~tt + , 
K~-k + tt°, and if~7r + 7r~7r + channels, where the mea- 
sured invariant mass is required to be within 20, 35, and 
20 MeV/c 2 , respectively, of the nominal D° mass [12], 
corresponding to 3.0, 2.5 and 3.0 a. A vertex fit is per- 
formed with the mass constrained to the nominal value. 
The 7T° candidate is formed from two photon candidates 
with invariant mass between 115 and 150 MeV/c 2 . 

For the measurement of the ratio r, the D° is combined 
with a low momentum ir to form a D* candidate, with its 
vertex constrained to the interaction point (beam spot). 
Candidates with mass difference mjja^ — m^n between 
144 and 147 MeV/c 2 are retained. A charged track, as- 
sumed to have the pion mass, is combined with the D* 
to form a B° candidate. The \ 2 probabilities for both 
the D* and B° vertex fits are required to be greater than 
0.1%. To reject jet-like continuum background, the nor- 
malized Fox- Wolfram second moment R2 [13], computed 
with charged tracks and neutral clusters, is required to 
be less than 0.5, and | cos^tI less than 0.85 where 9t is 
the thrust angle between the B° candidate and the rest 
of the event in the e + e~ center-of-mass (CM) frame. 

For B° -> D a K+7r- and D°K+tt- measurements, 
the B° candidate is formed by combining a D° can- 
didate with oppositely charged pion and kaon candi- 
dates. We select candidates outside the D*~K + region 
(142.5 < m D o n - m D a < 148.5 MeV/c 2 , a 6cr window). 
The measured D° invariant mass must be within 12, 28, 
and 8.5 MeV/c 2 of the nominal D° mass for Kir, Kirn , 
and Kttttt: modes, respectively. Candidates are rejected 
if the D° — > Kmr a decay probability, computed with the 
Dalitz parameters measured in Ref. [14], is less than 6% 
of the maximum value. The \ 2 probability of the D° 
(B°) vertex fit is required to be greater than 0.5% (2%). 
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All charged tracks are required to have at least 12 hits 
in the drift chamber and transverse momentum greater 
than 100 MeV/c. Both kaon candidates are required to 
be consistent with the kaon hypothesis. Prompt pion 
candidates consistent with the kaon hypothesis are re- 
jected. 

To further reduce the continuum background, | cos9g\ 
must be less than 0.9, where 9* B is the polar angle of 
the B° candidate in the CM frame. A Fisher discrimi- 
nant T is formed based on i?2 , cos 6t , 9* B , and two mo- 
ments Lq and L2, where Li — YljPjl cos 9* \ l , summed 
over the remaining particles j in the event, where 9* and 
p* are the angle with respect to the B° thrust and the 
momentum in the CM frame. Different cuts on T are 
applied for each mode to optimize the signal significance 
based on simulated event samples. Candidates used in 
the subsequent fits have beam-energy substituted mass 
toes = - O*) 2 > 5.2 GeV/c 2 and energy dif- 

ference \AE\ = \E* - yG/2\ < 150 MeV, where E* and 
p* are the energy and momentum of the B° candidate 
and ^/s is the total energy in the CM frame. 

We study five samples separately: (a) B° — > D°K + ir~ 
excluding the D*~ K+ contribution, (b) B° D°K + n~ , 
(c) B° -> D°K*{892)°, (d) B° -> D^(2460)-K+, and 
(e) B° — > D*~h + where h + is a pion or kaon. Samples 
(c) and (d) are subsets of (a), where the resonances are 
selected within 1.5 times their full widths [12]. 

For samples (a)-(d), a two-dimensional (toes, AE) 
unbinned-maximum-likelihood fit is used to determine 
the signal yields. The signal component is the product 
of a Gaussian in toes centered at the B° mass and a 
Crystal Ball lineshape [15] in AE centered near zero. 
The combinatorial background component is modeled 
with an Argus threshold function [16] in toes an d a 
second-order polynomial in AE. Two background com- 
ponents peak in toes: peaking background A describes 
the B° — > D**~ir + contribution, which also peaks in AE 
but the peak is shifted by about +50 MeV because the 
pion is misidentified as a kaon; peaking background B 
uses a second-order polynomial in AE to accommodate 
events such as D^K^tt, and D^*'p, where one or more 
pions or photons are missed in the reconstruction and/or 
a pion is misidentified as a kaon. The probability density 
function (PDF) is the sum of the signal and three back- 
ground components. A large B° — ► D*~tt + data control 
sample is used to determine the signal shape in both AE 
and toeSi an d the peaking background A in AE, where 
we assign the kaon mass to the pion candidate. We use 
the same parameters for signal and peaking backgrounds 
in Toes since they are consistent in simulation. The AE 
distributions and yields for the four components in the 
signal region are shown in Fig. 1 and Table I, respectively. 

The signal yield for B° — > D°K + ir~ is corrected for 
variations in signal efficiency across the DKtt Dalitz plot. 
Each event k with variables qu = (toes,A:, AEk) is as- 




-0.1 0.1 ° -0.1 0.1 

AE (GeV) AE (GeV) 



FIG. 1: AE distributions and PDF projections with r/iES > 
5.27 GeV/c 2 for (a) 5° -> WK + -k~ excluding D*~ K+ can- 
didates, (b) B° D°K + n-, (c) B° -> D°K* (892)° and (d) 
B° -> D5(2460)" K + , for the three D° modes combined. Cir- 
cles with error bars are data points. Four curves from top 
to bottom represent: the total PDF (solid), total background 
(dashed) , combinatorial background plus peaking background 
B described in the text (dot-dashed) and combinatorial back- 
ground only (dotted). In (a)-(c), the middle two curves over- 
lap because the peaking background A is negligible. 

signed a signal weight [17] 

Ei=i^ig,i p i(9fc) 
= -=3 — at p~7-T~\ » 

T,j=i N j p j(<ik) 

calculated from the four PDF components Pj , their yields 
Nj from the fit, and the covariance matrix elements V s i gj j 
between A^ig and Nj. The efficiency-corrected signal 
yield is then iy s ig(<f/c)/£fc, where the efficiency is 
estimated from the simulated events in the vicinity of 
each data point in the Dalitz plot. 

Figure 2 shows the signal weight distribution as a func- 
tion of mx+TT- an d m D°-K-- The peaks near TO^-.( 8 g2)o 
and Too* (2460)- are clearly visible. We use the (toes, AE) 
fit results and signal efficiencies estimated from simulated 
B° -> Z)°A'*(892)° and B° -> £> 2 , (2460)-AT+ samples 
to compute corresponding branching fractions. For the 
B° — > D°K + it~ mode, we assume a flat distribution on 
the Dalitz plot when determining the signal efficiency. 

For modes in which we do not observe a significant 
signal, the 90% confidence level (C.L.) branching frac- 
tion upper limit is determined by integrating the prod- 
uct of the PDFs for the three D° modes as a function 
of branching fraction from to 2?ul so that J Q UL CdB = 
0.9 J °° CdB, where C is the likelihood function. 

To measure r, we select events with toes > 
5.27 GeV/c 2 from sample (e). A two-dimensional PDF 
of AE and 9c (the reconstructed Cherenkov-light an- 
gle of the prompt track) is used to separate D*K from 
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TABLE I: The yields of signal, combinatorial (comb.) and peaking (peak A, peak B) background PDFs of the samples (a)- 
(d) described in the text; values and errors are rescaled to represent the yields in the signal region (tties > 5.27 GeV/c 2 , 
|Ai?| < 40 MeV). The bottom row shows the branching fractions with statistical errors. 



D° mode 


(a) B° -> D°K 


+ 7T" 


(b)B 


1 — > D°K 




(c)B 


-» D°JT (892)° 


(d)B 


3 -> Da (2460)^+ 


Kn Knn° 


Knnn 


Kn 


Knn° 


Knnn 


Kn 


Knn° Knnn 


Kn 
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Signal 


101 ± 17 58 ±20 


69 ± 19 


-17 ± 13 


34 ± 24 


8± 22 


35 ±7 


21 ±7 31 ±7 


15±6 


15±6 16±5 


Comb. 


229 ±4 500 ± 5 


528 ±5 


608 ±5 


918 ±6 


989 ±6 


17 ± 1 


29 ± 1 30 ± 1 


16 ±1 


16 ±1 22 ±1 


Peak A 


5±6 0±1 


0± 2 


0±0 


0±0 


0±0 


0±0 


0±0 0±0 


2± 2 


5±2 2±1 


Peak B 


45 ± 9 76 ± 12 


42 ± 10 


50 ±11 


54 ± 14 


45 ± 13 


6±3 


10 ±3 3±3 


2±3 


7±3 0±1 


B (1(T 6 ) 


88 ± 15 
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FIG. 2: The signal weight distribution as a function of 
m K + 7r - and m^jo^- ■ The shaded histograms include only 
events with (a) \m B o^- - 2460 MeV/c 2 | < 75 MeV/c 2 , and 
(b) \m K+n - - 896 MeV/c 2 | < 150 MeV/c 2 . 



D*tt decays. Tracks with an estimated 8c uncertainty 
<Jc > 4 mrad or n 7 . s /^/n 7 . s ± n 7i f, < 3 are removed, 
where n 7)S and n 7! b are the numbers of signal and back- 
ground photons determined from a likelihood fit to the 
ring of Cherenkov photons associated with the track [11]. 
Finally events are rejected if 8c is smaller than the pre- 
dicted Cherenkov angle for kaons by more than 4ctc7, in 
order to remove particles heavier than kaon. 

The AE signal peak PDF is a Crystal Ball lineshape 
and the background is a linear function plus a Gaus- 
sian peaked near —150 MeV to accommodate background 
events such as D*p and D**ir where a soft ir is missed 
in the reconstruction. The distribution of {6c — S c )/crc is 
modeled by Gaussian functions. For the pion component, 
we use three Gaussian functions centered near zero. For 
the kaon component, a single Gaussian function centered 
near [0q — 8 c )/ac is sufficient, where 8^ and 8 C are 
the expected Cherenkov angle for kaon and pion, respec- 
tively, based on the measured momentum. Most of the 
parameters are obtained from a fit to the pion or kaon 
tracks in a large cc D* X -> D°irX, D° -> R-tt+ 
data control sample, except the total width of the dis- 
tribution, which is free in the final fit to accommodate a 
small difference in width due to differences in momentum 
spectra between signal and control samples. 

Figure 3 shows the AE and (8c — 6 c )/ac distribu- 
tions and PDF projections for B° — > D*~h + (h = tt 



FIG. 3: (a) AE and (b) Cherenkov angle (9 c -0%)/ac distri- 
butions for D*~h + candidates and PDF projections. Circles 
with error bars are data points. Shaded distribution is com- 
binatorial background, the dotted curve adds the D*n con- 
tribution, and the solid curve is the full PDF. The dashed 
curve represents the D* K contribution only. AE for D*n is 
centered near zero, while for D* K it is shifted to lower values 
because the prompt track is assumed to be a pion. 

or K) candidates. We find 13400 signal events, of which 
/ = (6.80±0.28)% are D*K events, and 4850 background 
events in the sample. The ratio r = //(l — /) is cor- 
rected by the signal efficiency ratio r e — Ed*k / £d*-k — 
(94.0 ± 2.3)% obtained from simulation. This ratio is 
smaller than unity because 8c for kaons is smaller (re- 
sulting in fewer Cherenkov photons) and more kaons than 
pions decay in flight within the tracking volume. The un- 
certainty on r e includes simulation statistics and system- 
atic uncertainties due to the two aforementioned effects. 

For samples (a)-(d), the systematic uncertainties on 
the signal efficiency are studied with large r lepton de- 
cay samples (for track reconstruction efficiency) and com- 
parisons between signal simulation and B° — > D*~ir + 
data control sample. The fractional uncertainty, com- 
mon to all four samples, on signal efficiency is 5% in- 
cluding the uncertainties on the number of BB events 
and the D° branching fractions. For the B° — > D°K + ir~ 
mode, the uncertainty of efficiency variation on the Dalitz 
plot contributes an additional systematic error of 8%. 
In addition, we vary the control sample shapes in each 
fit by one standard error and sum the changes in sig- 
nal yield in quadrature. The total signal yield variations 
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are 8, 2.0, 3.4, and 2.6 events for D°K+ir-, D°K+ir-, 
D°K* (892)°, and D%(2460)-K+, respectively. For the 
B" -» D°K* (892)° and D?j(2460)-if + measurements, 
we consider possible contamination from each other and 
from the non-resonance contribution. Using the signal 
yields for B° -» D°X*(892)° and D^(24m)~K + , and 
the cross-feed efficiencies determined from simulation, we 
find that six events in each of these two B° modes could 
be attributed to the other mode and to non-resonance 
contributions. This contributes a 6% uncertainty for 
B° -> D°K*(892)° and 11% for B° -> D^{2AQQ)~K+. 
The uncertainty due to the full width of the Z^(2460) _ 
and K* (892)° resonances is 8% for B° -» D*(2460)-K+ 
and less than 1% for B° -> D°ir*(892) . 

The largest systematic uncertainties cancel in the 
branching ratio measurement (sample (e)). The remain- 
ing systematic errors are from PDF shapes, control sam- 
ple distributions and contaminations (1.9%), residual un- 
certainties in the signal efficiency ratio (2.4%), and po- 
tential fit bias (2.1%). The last item has been evaluated 
with simulation samples including background. 

In conclusion, we have measured the branching fraction 
for the B° -> D°K+7r- decay excluding D*~ K+ , 

B(B° -» D°K+ir-) = (88 ± 15 ± 9) x 10~ 6 , 

as well as its two significant resonances, 

B(B° -> D°K*(892)°) ■ B(K*(892)° -» K+n') 
= (38 ± 6 ± 4) x 10~ 6 , and 
B(B° -» D*(2460)~ K+) ■ £(D*(2460)- -» D°tO 
= (18.3 ±4.0 ±3.1) x 10~ 6 . 

The signal significances are 8.7, 8.3 and 5.0 standard de- 
viations, respectively, determined from the change in the 
likelihood between the best fit and a fit with the signal 
yield fixed to zero (the first case) or the possible cross 
feed from other sources (six events for the latter two 
cases). From a fit excluding the observed resonances, 
assuming flat distriubtion on the Dalitz plot, we find 
B(B° -» D a K+ir-) = (26 ± 8 ± 4) x 10~ 6 , whose sig- 
nal significance is 3.1a and 90% confidence level upper 
limit is 37 x 10~ 6 . We do not observe a significant sig- 
nal for the CKM-suppressed B° -> D°K+ir- mode. The 
90% confidence level upper limit is B(B° -> D q K + -k~) < 
19 x 10~ 6 . The event yields in this channel are lower than 
anticipated [5] indicating that a significantly larger data 
sample is required to constrain 7 through this method. 

The ratio of branching fractions for B° — > D*~K + to 
B° — > D*~n + is measured to be 

r = (7.76 ±0.34 ±0.29)% , 

a nearly four-fold improvement compared to the pre- 
vious result [8]. This ratio is consistent with 
(/W/7r) 2 tan 2 6>cab - °- 072 I 18 ] j expected at tree level 



if factorization and flavor-SU(3) symmetry hold, where 
#Cab is the Cabibbo angle and fx and f„ are the decay 
constants of the kaon and pion, respectively. 
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